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Abstract 

We present analysis of optical photometric and polarimetric observations of contact binary system FO Hydra (FO 
Hya). The computed period of the system is 0.469556±0. 000003 days. An O-C curve analysis indicates an increase 
in its period by 5.77 x 10~ 8 day yr _1 . The photometric light curves are analyzed using Wilson-Devinney code (WD). 
The present analysis shows that FO Hya is a B-subtype of W UMa type contact binary. The radii and mass of primary 
and secondary components are found, respectively, to be Ri = 1.62±0.03 R Q and R2 = 0.91±0.02 R , and Mi = 
1.31+0.07 M and M2 = 0.3 1±0. 1 1 M . The light curve shape shows small asymmetries around the primary and 
secondary maxima. This may be due to the presence of dark spots on the components. The polarimetric observations 
yield average values of its polarization to be 0.18±0.03, 0.15+0.03, 0. 17+0.02 and 0.15+0.02 per cent in B, V, R and 
I bands, respectively. These values are appreciably lower than the typical polarization of W UMa type binaries. We 
have discussed the possible sources of the observed polarization but in order to arrive at definitive answer to this more 
phase locked observations are needed. 
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1. Introduction 

Eclipsing binary stars of W Ursae Majoris (W UMa)- 
type variables; referred to, in short, EW, have their 
light curves with strongly curved maxima and minima 
that are nearly equal in height/ depth. These stars are 
divided into two sub-classes: A-type and W-type. In 
the A-type systems the larger component has the higher 
temperature whereas in the W-type systems the smaller 
component has the higher temperature. However, the 
temperature difference in the primary and the secondary 
is less than 1000K. Observationally, it has been found 
that the A-type systems tend to have low mass ratios (q 
< 0.3) and spectral type from A to F. W-type systems 
usually have mass ratios q > 0.3 and spectral types 
of G or K. Besides, there is also another subclass 
of contact binaries having temperature difference of 
1000 K or larger between their components (Rucinski, 
1986; Eggleton, 1996). This temperature difference 
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is considerably larger than the theoretical limit of 
about 1000 K (e.g. Kaluzny, 1986; Hilditch and King, 
1988). Lucy and Wilson (1979) introduced the class 
of B-type systems which are systems in geometrical 
contact, but not in thermal contact and therefore there 
are large surface temperature differences between the 
components. The B-type systems are sometimes re- 
ferred to as poor thermal contact systems (Rucinski and 
Duerbeck, 1997). Further, the light curves of W UMa 
binaries usually show differences in brightness of their 
maxima; this asymmetry is called the O'Connell effect 
(O'Connell, 1951). The O'Connell effect has been of 
particular interest in understanding the light variations 
of W UMa binaries. It also has been reported that the 
shape of the light curves and the depth of eclipses vary 
with time for some systems. Thus, the asymmetric light 
curves of W UMa type systems indicate the possible 
presence of starspots on one or both stars of the system. 

The eclipsing binary FO Hydra (BD-18 2822) 
has been discovered to be a variable by Hoffmeister 
(1936). First quantitative observations were made by 
Tsessevich (1954) who reported a period of 1.159 
days. The later studies of FO Hya indicate a slightly 
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asymmetric light curve having strong curvature in both 
maxima and minima of unequal depths (Binnendijk, 
1960; Kopal, 1978; Candy and Candy, 1997), which 
was explained by the presence of spot on a star of 
the binary system. Candy and Candy (1997) have 
suggested that the slight and irregular variation in the 
orbital period of FO Hya is presumably caused by mass 
transfer between the components. They determine the 
orbital period to be 0.4695571 days. Recently, the 
light and radial velocity curve analysis was performed 
by Siwak et al. (2010). They find the contact nature 
of FO Hya with temperature difference of ~ 2333 K 
between components. However, their light curves have 
been incomplete at the first quadrature and because of 
this, they were unable to get O'Connel effect and other 
parameters exactly. Therefore, we have undertaken 
an investigation of this system in order to get a better 
parameters estimation. 

In this paper, we report broad band B, V, R and I light 
curves analysis of the star FO Hya using our photomet- 
ric observations. We describe the observations and data 
reduction procedures in §2. In §3 we present the pe- 
riod analysis. The light curves and estimation of basic 
parameters are subject matter of §4 and §5 . In §6, we 
present the analysis of the observations using Wilson- 
Devinney light curve modeling technique (WD code). 
We discuss the sources of polarization in §7. The sum- 
mary of the results is given in the last section. 

2. Observations and Data Reduction 

2.1. Optical Photometry 

The photometric observations during 11 nights be- 
tween April 09, 2009 to March 04, 2010 have been 
carried out using broad band B, V, R and I filters at 
Aryabhatta Research Institute of observational sciencES 
(ARIES), Nainital with a 2kx2k CCD camera mounted 
on cassegrain focus of the 104-cm Sampurnanand tele- 
scope (ST). The exposure times are ranging from 10 to 
90 seconds depending upon the individual system and 
the observing conditions. The CCD system consists of 
24x24 fi 2 size pixel. The gain and readout noise of CCD 
are 10 e-/ADU and 5.3 e-, respectively. The CCD cov- 
ers a field of view ~13xl3 arcminute 2 . Several bias and 
twilight flat frames were also taken during each observ- 
ing run. Bias subtraction, flat fielding^and aperture pho- 
tometry were performed using IRAI 



We have chosen 2mass 09594720-1908422 and 2mass 
09594342-1905435, as the comparison and check stars, 
respectively. The image of the field depicting the com- 
parison and the check stars with FO Hya at the center; 
taken from ESO online Digitized Sky Survey, is shown 
in Figure 1. The co-ordinates of the target, compari- 
son and the check stars are listed in Table Q] Differ- 
ential photometry, in the sense of variable minus com- 
parison star, have been performed and the comparison- 
check light curve in V band along with the folded V- 
band light curve of FO Hya is shown in Fig. 2. The 
ephemeris used is as given by equation (1) of the §3. 
From the light curve, it is clear that the comparison and 
check stars are constant thought the observations. The 
nightly mean of standard deviation (<x) of different mea- 
sures of comparison minus check stars in B, V, R and I 
band are 0.019, 0.013, 0.012 and 0.01 1, respectively. 
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'IRAF is distributed by National Optical Astronomy Observato- 
ries, USA; |http://iraf.net| 



Figure 1: A ~ 10' xlO' image of the field with FO Hya at center is 
taken from ESO online Digitized Sky Survey. The target, comparison 
and the check stars are marked. 



2.2. Optical Polarimetry 

The broad-band B, V, R and I polarimetric obser- 
vations of FO Hya have been obtained on April 07, 
23, 2010 and May 20, 2010 using ARIES Imaging 
Polarimeter (AIMPOL) mounted on the cassegrain fo- 
cus of the 104-cm Sampurnanand Telescope of ARIES, 
Nainital. The detail descriptions about the AIMPOL, 
data reduction and calculations of polarization, position 
angle are given in Rautela et al. (2004). For calibration 
of polarization angle zero-point, we observed highly po- 
larized standard stars, and the results are given in Table 
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Table 1: Coordinates and magnitudes of the program (FO Hya), comparison (2Mass 09594720-1908422) and check (2Mass 09594342-1905435) 
stars. 



Star 


RA (J2000) 


Dec (J2000) 


B mag 


V mag 


J mag 


FOHya 


09 59 43.29 


-19 07 56.22 


11.31 


11.05 


10.36 


2Mass 09594720-1908422 


09 59 47.21 


-19 08 42.27 






11.41 


2Mass 09594342-1905435 


09 59 43.43 


-19 05 43.54 






11.99 



[2] Both the program and the standard stars were ob- 
served during the same night. The obtained values of 
polarization and position angles for standard polarized 
stars are in good agreement with Schimidt et al. (1992). 
To estimate the value of instrumental polarization, un- 
polarized standard stars were observed. These measure- 
ments show that the instrumental polarization is below 
0.1 per cent in all pass bands (see Table 0. In fact, 
the instrumental polarization of ST has been monitored 
since 2004 within other projects as well (e.g. Medhi et 
al., 2007, 2008; Pandey et al., 2009). These measure- 
ments demonstrated that it is invariable in all B, V and 
R pass bands. The instrumental polarization was then 
applied to all measurements. 

Table 2: Observed polarized and unpolarized standard stars 
Polarized Standard 



Schmidt et. al(1992) 



Present work 



Filter 


P (per cent) 


en 


P (per cent) 


en 






HD 154445 






B 


3.445 ± 0.047 


88.88 ±0.39 


3.836 ±0.211 


88.88 ± 1.54 


V 


3.780 ± 0.062 


88.79 ± 0.47 


3.708 ± 0.098 


88.79 ± 0.78 


R 


3.683 ± 0.072 


88.91 ±0.56 


3.652 ± 0.029 


88.91 ±0.22 


I 


3.246 ± 0.078 


89.91 ±0.69 


3.353 ± 0.006 


89.91 ±0.05 






HD 155197 






B 


4.1 12 ±0.047 


103.06 ±0.33 


4. 109 ±0.039 


103.06 ±0.27 


V 


4.320 ±0.023 


102.84 ±0.15 


4.283 ± 0.006 


102.84 ±0.04 


R 


4.274 ± 0.027 


102.88 ±0.18 


4.285 ±0.037 


102.88 ±0.24 


I 


3.906 ± 0.041 


103.18 ±0.30 


3.679 ± 0.060 


103.18 ±0.45 






HD 154445 






B 


3.445 ± 0.047 


88.88 ±0.39 


3.640 ± 0.084 


88.88 ± 0.65 






HD 161056 






B 


3.799 ± 0.055 


66.56 ± 0.42 


3.751 ±0.055 


66.56 ±0.41 




Unpolarized Standard (/? Uma) 








q(per cent) 


u(per cent) 






B 


-0.114 


-0.091 






V 


- 0.264 


-0.133 






R 


- 0.001 


- 0.070 






I 


0.005 


- 0.079 







3. Period analysis 

There are number of techniques available for study- 
ing the periodicity in a time series data. We have derived 
the period of FO Hya using CLEAN algorithm (Roberts 
et al., 1987) method. The power spectrum is shown in 
Fig. 3 wherein the peak frequency exits at ~ 4.2592 
(day -1 ). Eclipsing binary light curves can be repre- 
sented by two sine waves and therefore multiplying the 



period (=l/frerquency) so obtained by 2 gives appropri- 
ate orbital period of the binary system. We determine 
the orbital period of FO Hya to be 0.46957±0.00004 
day, which is in agreement with the period derived by 
Candy and Candy (1997) and Siwak et al. (2010). 
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Figure 2: Folded V-band light curve of FO Hya along with the 
Comparison-check folded light curve in V band. The ephemeris is 
given in equation 1 . 



The timings of light minima are derived using Kwee 
and van Woerden's (1956) method. We derive the tim- 
ings of the 10 minima, using our quasi-simultaneous B, 
V , R and I bands observations and of the 25 minima, 
using All Sky Automated Survey (ASAS; Pojmanski, 
2002) data. We also use the world-wide database of 
minima timings of Paschke and Brat (20062 and obtain 
a total of 67 minima timings over a span of 65 years. 
The results are given in Table [3] 

We determine the ephemeris as given by the follow- 
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4 5 6 
Frequency (c/d) 

Figure 3: Power spectra of FO Hya obtained from CLEAN algorithm. 
The most dominant peak in the spectra occurs to ~4.259 c/d. 



xlO" 8 day yr -1 , which is equivalent to 7.42 xKT 11 
day cycle 1 ; meaning there by a period increase of 
0.50 s century" 1 (see Qian and Ma (2001) for the 
details of the calculation). The change in the period is 
not unusual among W UMa type stars. The observed 
period increase/decrease rates of these W Uma type 
binaries are in the range from ±10~ 9 to + 10~ 6 day yr -1 
(Dryomova and Svechnikov, 2006). Several W UMa 
systems show the increasing trend in their orbital period 
e.g. YY Eri (Kim et al., 1997), V839 Oph (Akalin and 
Derman, 1997), AB And (Kalimeris et al., 1994), DK 
Cyg (Awadalla, 1994), V401 Cyg (Herczeg, 1993), 
V566 Oph (Maddox and Bookmyer, 1981). Examples 
of W UMa systems showing decreasing trend in their 
orbital period are TZ Boo, Y Sex (Qian and Liu, 2000), 
V502 Oph (Derman and Demircan, 1992), U Peg (Zhai 
et al., 1984), VW Cep (Kaszas et al., 1998) and V781 
Tau (Liu and Yang, 2000a). 



ing equation 

Min.I(HJD) = 2455259.32959(64) + d .469556(3) xE 

(1) 

The numbers given in parentheses represent the proba- 
ble errors and are expressed in terms of the last quoted 
digits. Now we compute the O-C residuals, and plot 
its variation against the epoch number in Fig. 4 . As 
seen from the figure the general trend of the (O-C) i vs 
epoch diagram indicates a roughly parabolic distribu- 
tion indicating a long time change in the orbital period. 
We, therefore, attempt second order polynomial to fit 
the data and represent it by a dashed line in the top panel 
of the Fig. 4. Here we also represent the residual from 
polynomial fit (O-Q2 in the middle panel of the Fig. 
4. It is noticed that there is still another change, which 
seemingly may be sinusoidal. The continuous curve in 
the middle panel of the Fig. 4 shows the best sinusoidal 
fit. We also plot the residuals of the sinusoidal fit, shown 
in the bottom panel of the figure, indicating no signifi- 
cant variation. The combined quadratic and sinusoidal 
fit for the observed data is obtained as: 

O-C = 0.0045(+0.0013) 
+(1.01 +0.18) x 10~ 6 x£ 
+(3.71 + 0.39) x 10~ n xE 2 (2) 
+0.0041(±0.0018)sw[(1.91 + 0.18) x 10~ 4 x E 

+3.93(±0.31)]. 

The coefficient of the square term yields the rate of 
change of the period. We determine (dP/dt) to be 5.77 



As the two components of FO Hya are in a state of 
over-contact, we take the most likely cause of the pe- 
riod change as the transfer of mass between the compo- 
nents. The amount of mass transfer, dm is related to the 
change in period, dP, for a system of total mass, M via 
(Coughlin et al., 2008; Yang and Liu, 2003). 

dm Mq dP 

~dt ~ 3P(1 - q 2 ) ~dt 

By employing M = 1.62+0.09 M G and q = 0.238; the 
values obtained in section |6l mass transfer rate of FO 
Hya between the components turns out to be 1.67 x 
lO^Moyr- 1 . 

The sinusoidal term of the (O-C) oscillation in equa- 
tion d2]i reveals a periodic change with semi amplitude, 
A, of 0.0041 days and a period, P3 of 42.30 years. The 
amplitude of period oscillation, AP may now be de- 
termined using the following relation due to Rovithis- 
Livaniou et al. (2000). 

AP = V2[l - cos(2ttP/P 3 )] x A (4) 

The amplitude of the period oscillation is obtained to 
beAP = 0.78 xl0~ 6 day leading to AP/ P = 1.66 xl0~ 6 . 

The cyclic variation may be interpreted as resulting 
either from the magnetic activity of one or both com- 
ponents (Applegate, 1992), or from the light-time ef- 
fect (LITE) through the presence of a tertiary compan- 
ion (Irwin, 1952, 1959). The required variation of the 
quadruple momentum AQ so as to reproduce this cyclic 
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Table 3: Times of light minima of FO Hya. 



HJD 

(2400000+ ) 


rr 

Type 1 


Method 2 


Ref/ 


HJD 

(2400000+) 




Type 1 


Method 2 


Ref. J 


31216.17000 


P 


vis 


(1) 


53123.54854 


s 


ccd 


(7) 


31224.21000 


P 


vis 


(1) 


53129.65254 


s 


ccd 


(7) 


31231.23000 


P 


vis 


(1) 


53133.64468 


p 


ccd 


(7) 


44618.23440 


P 


pe 


(2) 


53393.78252 


p 


ccd 


(7) 


44620.11110 


P 


pe 


(2) 


53400.82116 


p 


ccd 


(7) 


44632.09000 


s 


pe 


(2) 


53417.72362 


p 


ccd 


(7) 


44632.32320 


P 


pe 


(2) 


53436.74512 


s 


ccd 


(7) 


44634.20050 


P 


pe 


(2) 


53499.42900 


p 


V 


(8) 


44638.19000 


s 


pe 


(2) 


53500.60707 


s 


ccd 


(7) 


44640.30510 


P 


pe 


(2) 


53504.59680 


p 


ccd 


(7) 


44642.18270 


p 


pe 


(2) 


53556.48409 


s 


ccd 


(7) 


44649.22710 


p 


pe 


(2) 


53673.86840 


s 


ccd 


(7) 


44656.26980 


p 


pe 


(2) 


53720.82771 


s 


ccd 


(7) 


44659.08900 


p 


pe 


(2) 


53724.81814 


p 


ccd 


(7) 


46137.24970 


p 


pe 


(2) 


53792.67499 


s 


ccd 


(7) 


46138.19050 


p 


pe 


(2) 


53836.57843 


p 


ccd 


(7) 


46144.07000 


s 


pe 


(2) 


53867.56291 


p 


ccd 


(7) 


46144.29840 


p 


pe 


(2) 


53871.56083 


s 


ccd 


(7) 


47970.18000 


s 


pe 


(2) 


53903.48665 


s 


ccd 


(7) 


47971.11000 


s 


pe 


(2) 


54115.26260 


s 


Ic 


(9) 


47974.15940 


p 


pe 


(2) 


54157.75581 


p 


ccd 


(7) 


48300.27000 


s 


pe 


(2) 


54464.84950 


p 


ccd 


(7) 


48301.21000 


s 


pe 


(2) 


54505.69046 


p 


ccd 


(7) 


48302.14000 


s 


pe 


(2) 


54932.28417 


s 


ccd 


(10) 


48303.32030 


p 


pe 


(2) 


54933.22186 


s 


ccd 


(10) 


50189.00300 


p 


vis 


(3) 


54934.16650 


s 


ccd 


(10) 


50189.00300 


p 


V 


(4) 


55254.16519 


p 


ccd 


(10) 


51286.87600 


p 


ccd 


(5) 


55254.39782 


s 


ccd 


(10) 


51905.28800 


p 


V 


(5) 


55255.33811 


s 


ccd 


(10) 


52347.14330 


p 


ccd 


(6) 


55257.22290 


s 


ccd 


(10) 


52960.85457 


p 


ccd 


(7) 


55258.39126 


p 


ccd 


(10) 


53079.64647 


p 


ccd 


(7) 


55258.16047 


s 


ccd 


(10) 


53101.71414 


p 


ccd 


(7) 


55259.32907 


p 


ccd 


(10) 


53109.69754 


p 


ccd 


(7) 











Type: p = primary minimum lime, s = secondary minimum lime. 
2 Method: vis - visual observation, pe = photoelectric data, V = Johnson V filter , ccd = charge- coupled device, Ic = Cousins I filler. 

References: (l)Zessewiisch, V. P. (1954); (2) Candy and Candy (1997); (3) Nagai, Kazuhiro; (4) Nagai, Kazuo; (5) Paschke, Anton; (6) Kiyota, Seichiiro; (7) ASAS; Pojmanski (2002) (8) Sobotka, Petr (2007); (9)Nakajima, Kazuhir; 
(10) Present Study. 
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change, may be calculated using the following equation 
due to Lanza and Rodono (2002). 



P Ma 1 



(5) 



Using equation © we get AQ as 2.08 xlO 49 and 0.49 
xlO 49 g cm 2 for primary and secondary components, 
respectively. The AQ values are smaller than its typical 
value of 10 s 1 - 10 52 g cm for close binaries (Lanza 
and Rodono, 1999), suggesting, that the magnetic 
mechanism may not suffice significantly to explain the 
cyclical period variation. Hence it turns out that the 
cyclical period change is more plausibly be interpreted 
as due to the presence of a third body. 

The projected orbital semi- major axis, a i2sin//, of the 
orbit of the third body is related to the semi-amplitude 
of the O-C oscillation via ansiai/ — Ax c (Rovithis- 
Livaniou et al., 2000), where if is the inclination of the 
orbit of the third component and c is the speed of light. 
We obtain a^sim'/ to be 0.71 AU. Now, it is straight- 
forward to obtain the mass function of the third body, 
f(M 3 ) using the following well-known equation. 
Here a represents the separation of the components. Let 
us recall that FO Hya is a spectroscopic binary and 
hence its radial velocity curve analysis, which has been 
made earlier by Siwak et al. (2010), may also be used 
to determine 'a'. The well known radial velocity ampli- 
tude parameters Ki and K2 are related to the respective 
ai and a2 via 



Ki(2) - 



2n ai(2)Sim 



(6) 



Making use of the values obtained by Siwak et al. 
(2010); Ki = (62.48±0.97) km/s, K 2 = (253.20±4.50) 
km/s and of the period and inclination as obtained by 
us we get a! = (0.59+0.02) R Q , a 2 = (2.40±0.06) R 
leading to a = (2.99+0.06) R Q and ai/a 2 = M2/M1 = 
q = 0.238. Now, using Kepler's third law it is straight 
forward to determine M and thereby to obtain Mi 
1.31-0.07 XI and M 2 =0.31 ±0.11 M Q . 



An 

f(M 3 ) = — 2 x(fli 2 sin//) 3 



(7) 



We obtain the mass function of the third body to be 
f(M 3 ) = 0.0002 M . The mass function /(M 3 ) may 
also be expressed via 



f(M 3 ) 



(M 3 sin if) 3 
(Mi +M 2 + M 3 ) 2 



(8) 



where M\, M 2 , and M3 represent the masses of the 
eclipsing pair and the third companion, respectively and 
G is the gravitational constant. Now if one assumes the 
third body orbit to be coplanar to the orbit of the eclips- 
ing pair (i.e., it = 78.°05(+0.°43), the value of the low- 
est mass of the third body is computed to be 0.10 M . 



4. Light curves from present observations 

Fig. 5 shows the folded light curves of FO Hya in 
B, V, R and I bands. A rounded bottom primary (phase 
= 0.0) and almost fiat secondary (phase = 0.5) eclipses 
were seen which could be due to a total-eclipse config- 
uration of the system. The primary eclipse is obviously 
deeper than the secondary one. The difference between 
the primary and secondary minima is found in the range 
0.26 - 0.30 magnitudes in B, V, R, and I bands. There- 
fore, primary minimum is deeper and occurs when the 
hotter star passes behind the cooler one. On the other 
hand, shallower secondary minimum occurs when the 
cooler star passes behind the hotter one. The difference 
between primary maximum (phase = 0.25) and primary 
minimum is about 0.50, 0.48, 0.48 and 0.46 magnitudes 
in B, V, R and I bands, respectively. However, the differ- 
ence between secondary maximum (phase = 0.75) and 
secondary minimum ranges from 0.23 to 0.24 mag in B, 
V, R, and I bands, respectively. Similar differences be- 
tween maximum and minimum has also been reported 
earlier by Candy and Candy (1997) and Siwak et al. 
(2010). 

In W UMa-type systems, as both the stars are very 
close to each other, so, there is a continuous light varia- 
tion outside the eclipses. Moreover, the stars will expe- 
rience gravitational distortion and heating effects. These 
effects shows the asymmetric light curves of W UMa- 
type systems. The observed light curves of FO Hya 
are found to be asymmetric around the primary and sec- 
ondary maxima, with the secondary maximum brighter 
than the primary maximum. The magnitude difference 
between phase at 0.75, which is brighter, and phase at 
0.25 ranges from 0.039 to 0.041 mag in B, V, R, and I 
bands. This is typically explained by the presence of a 
spots on the components of the binary system. This in- 
dicates that light curves of FO Hya shows the O'Connell 
effect. The O'connel effect has also been noticed in 
many W UMa-type binaries e.g. FG Hya (Qian and 
Yang, 2005), CU Tau (Qian et al., 2005), AH Cnc (Qian 
et al., 2006), AD Cnc (Yang and Liu, 2002a), QX And 
(Qian et al., 2007), CE Leo (Yang and Liu, 2002b; Kang 
et al., 2004), BX Peg (Lee et al., 2004). 
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5. Estimation of basic parameters of FO Hya 

We determine the effective temperature of the pri- 
mary component employing the results due to Wang 
(1994) on contact binaries. First of all, we compute the 
color of the system using their equation (17), period- 
color relation, given as : 



(B- V)o = 0.077 -l.003logP 



(9) 



where orbital period P is in days. This relation yields the 
intrinsic color of the star system to be (B-V)o = 0.406. 
Now we use this value of (B-V)o in their Eq. (14) i. e. 



(B - V) = 



3.970 - logT ef f 
0.310 



(10) 



and get T e ff of the primary component to be 6982 K. 
On the other hand, the infrared color index (J-K), using 
2MASS catalogue (Cutri et al., 2003) turns out to be 
0.235. This value when used in the color-temperature 
calibration of Cox (2000) yields a temperature of 7000 
K and a spectral type F0V for the primary component 
of the binary system. The interstellar reddening is 
negligible in near infrared (Pribulla et al., 2009). 
Therefore, {J - K) color from 2MASS is reliable to 
derive the spectral type of FO Hya. 
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Figure 4: The top panel: (O-C)i diagram computed using equation 
(1). The general trend of the (O-C)i curve reveals a long-term period 
increase (dashed line). The middle panel: the (0-C)2 curve from the 
parabolic fit and its description by a sinusoidal equation (solid line). 
The bottom panel: residuals from equation (2). 



Rucinski and Duerbeck (1997) has obtained a calibra- 
tion of absolute magnitude (My) of W UMa-type star 
systems in terms of their (B - V)o color and period as 

M v = -AAAlogP + 3.02(B - V) + 0.12 (11) 

We obtain My of FO Hya to be 2.80 mag. Taking the 
visual magnitude (V) of FO Hya as 1 1 .05 and assuming 
the negligible reddening, the distance modulus of FO 
Hya is calculated to be 8.25 mag leading to a distance 
of 446.68 pc for the FO Hya. 

6. Light Curve Modeling and Photometric Solutions 

We have modeled B, V, R and I band light curves us- 
ing the WD code (Wilson and Devinney, 1971) imple- 
mented in PHOEBE^ (Prsa and Zwitter, 2005). It is a 
modified package of the widely used WD program for 
deriving the geometrical and physical parameters of the 
eclipsing binary stars. In the WD code, some of the pa- 
rameters need to be fixed and the convergence of the 



3 http://phoebe.fiz.uni-lj.si/ 



solutions are obtained using the methods of multiple- 
subsets (Wilson and Biermann, 1976). Since most of 
the W UMa binaries are the contact binaries, therefore, 
mode-3 of WD code is the most suitable for the light 
curve analysis. However, we start with mode-2 (for de- 
tached binary) and achieve the convergence of the so- 
lution and thereafter go to mode-3 for the final conver- 
gence. In the light curve analysis using the WD code, 
star 1 is the one eclipsed at primary minimum and star 
2, at the secondary minimum. This notation has been 
used throughout the modeling. 

The temperature of the star 1, T; = 7000K, the bolo- 
metric albedos Ai = A2 = 0.5 for convective envelopes 
(Rucinski, 1969) and gravity darkening coefficients, gi 
= g2 = 0.32 for convective envelopes (Lucy, 1967) are 
input parameters. These have been kept fixed to ob- 
tain exact solution. The bolometric (X\boi, y\boi, ^iboh 
yiboi) and monochromatic (xj, X2, yi, yi) limb darken- 
ing coefficients of the components are interpolated us- 
ing square root law from van Hamme (1993) tables. The 
adjustable parameters have been the temperature of the 
star 2 (T2), orbital inclination (i), the surface potentials 
of both components, Q.\ and Q2; Q.\ = Q.2 (for contact 
binaries), and monochromatic luminosity of star 1, Lj. 
During the photometric solution, a third object, with 
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Table 4: Photometric solutions obtained from the light curve modeling of the FO Hya. 



Parameters 


Symbol 


Photometric elements 


Mass Ratio iMi/Mii 


4 


238+0 006 


Orhital iiipliiiatioTi r°l 


j 


78 05+0 43 


I? £i 1 o ti \t'(^ lnminncitv nf ctni* 1 [R 1 

IVCltlLlVC lUllllllUalLy Ul iLd.1 1 I 1_> ) 


n . /CT , + T -.yi 




\ v ) 


IMAM ~ MU 


0.923 




IMAM + M^J 


Q21 


(T) 


rr ,/rr 1 + T itl 

IMAM ^ 


908 


IVlUllU^IllUlllallL. 1LI V JJ ^ 


v A lB — A 25/ 


1 Q? 




vyiB — 


fSQO 


1. * ) 


v A lV — A 2Vj 


0fi3 




vyiv — y2v7 


0.724 


fR") 


f Y 1 n Yl 


044 




\yiR — yiR) 


W. / JU 


W 


(x 1 i — Y-i 


0.105 




{V* i V-i 

Ly 1/ — yiij 


698 

U.U7U 


Rnlnmptrip 1H popffipipritG 


(x i ; / — Y-w A 


0.085 




vyiBo/ — yiboij 


0.638 


Third lisht (B) 

A. 1111 ^1 lliLllL V 1— x 1 




05901+0 00003 


V v ) 


I3V 


05169+0 00002 


(R) 


t3fi 


03367+0 00009 




l.. 


03454+0 00003 

U.uJI^JtIU.UUUUJ 


^lirfapp Potential*; 

OLllltl^t^ 1 VJLC^11L1H1l3 


v ail — ±L2) 


2 221 1 +0 0001 


J.1111C1 ^VJllLtl^L k>Lll itl^C 


(CI "1 


2.3247 




en a 


2.1742 


Fill-nnt faptnr 

J. Ill VJLIL laA^LVJl 


f 


0.68 


Kffective temneratures 

1 — - 11 V. V LA V w IV 1 1 1 l-> \~ 1 Ll L Ll 1 U'J 


Ti (K) 


7000 




L 2 \**-) 


^21 3 1+02 

JZ.1 J. J. _!_ V . 


v.l Ll 1 X IIL V Ll 1 L / V_l V V l3 \ UUlUlllL' L 1 Awl 


A. = At 

-^2 


0.50 


frravitv HarlrPTiiTiCT f^npffipipiit*; 

VJiav il y uui avI/1aaii£i v^v^l- i i il. i l^ 1 1 l.i 




0.32 


^viiphrnriipitv nnrmnptprs; 
oy 11^111 vjiii^il y i j tii tii i il Lti a 


Fi — F2 


1.0 


Rplativp RnHii nf pnmnnriprits; 

IVC^lULl V t JA.LIU11 V ) I ^VJlllJJVJllC^llLj 








ri 1 

1 1 pole 


498+0 001 






0.272+0.002 




^ I side 


547+0 002 




^2side 


0.288+0.002 


(back) 


*lback 


0.579+0.002 




*2back 


0.356+0.006 


Spot parameters 






Spot latitude 


[0(rad)] 


1.571 


Spot longitude 


[O(rad)] 


1.289 


Spot radius 


[r. s (rad)] 


1.078 


Spot temperature factor 


(T s /T.) 


0.890 


Absolute bolometric magnitudes 


M 6o/ i 


2.77 




M io;2 


3.03 


Effective gravity of the components 


log gi (cgs) 


4.14+0.07 




log g 2 (Cgs) 


4.01 + 1.25 


Semi-major axis 


[a or »(R.sim)] 


2.99+0.06 
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light I3, has also been taken as an adjustable parame- 
ter in order to get better fit of the light curves. The third 
light I3 is in the unit of total light as implemented in 
PHOEBE. However, determination of third light from 
the photometric light curve modeling sometimes affects 
the orbital inclination and amplitude of the light curve. 
But we can not ignore the third light I3 for our present 
light curve modeling, because from our O-C curve anal- 
ysis (see section [3]), we have found the evidence of a 
third light 1 3 . 

In order to derive reliable physical parameters of bi- 
nary stars from photometric light curve modeling one 
requires the spectroscopically determined mass ratio 
(Deb and Singh, 2011). The spectroscopic mass ratio 
q spec = 0.238 (Siwak et al., 2010) of FO Hya has been 
adopted during the modeling of light curves. We find 
that the theoretical light curves do not fit very well the 
observed ones; especially around the maxima. There- 
fore, the spot models are adopted to fit the observed 
light curves. We consider the possible spot models of 
hot or cool spots on the primary and secondary com- 
ponent, respectively. We then carry forth by iteration, 
allowing the spot latitude (ranging from 0° at the north 
pole to 180° at the south pole), the longitude (ranging 
from 0° to 360°, with 0° at the inner Lagrangian point, 
180° at the back end, and increasing in the direction of 
rotation), the angular radius (where 90° covers exactly 
half the star), and the temperature factor (the ratio of the 
spot temperature to the underlying surface temperature), 
as well as T2 to vary, until a satisfactory fit is found until 
any further corrections were less than the errors. At this 
point, having obtained both orbital and spot parameters, 
we solve it further for a final solution for folded un- 
binned light curve, allowing all previously mentioned 
parameters to vary, once again until any further correc- 
tions were less than the errors. It turns out that placing 
the spot on the secondary component is found to be bet- 
ter than the spot on primary. Therefore, we take this as 
our final (best fit) solution that corresponds to the min- 
imum sum of the square of the residuals value among 
them. The best fit theoretical light curves (solid lines) 
along with the observed light curve are shown in Fig. 5. 
The corresponding geometric configurations of FO Hya 
with a dark spot on the secondary component are plot- 
ted in Fig. 6 at phases (0) = 0.00, 0.25, 0.50 and 0.75, 
respectively. 

The spot solutions included in the WD code can 
be easily fitted to the whole light curve but poses a 
serious problem regarding uniqueness of the solution 
(Maceroni and van't Veer, 1993) unless other means of 
investigation such as Doppler Imaging techniques are 




0.0 0.2 0.4 0.6 0.8 1.0 

Phose (0) 



Figure 5: Phased light curves of the FO Hya in B, V, R and I broad 
bands, respectively. Open circles denote the observational data points. 
The continuous line is the synthetic light curves computed from the 
WD light curve modeling technique considering the case as presence 
of spot on the secondary component. 



applied (Maceroni et al., 1994). The WD's differential 
correction (dc) minimization program yields the values 
of the fitting parameters as well as the formal statistical 
errors associated with each of them. We perform the 
Monte Carlo parameter scan (heuristic scan) around 
the best solution making use of the PHOEBE's scripter 
capability (Bonanos, 2009) with a view to explore 
the values, errors and stability of the solutions. We 
run the WD's (dc) minimization program 1000 times; 
updating each time the input parameter values for 
the next iteration. We obtain the final values for the 
parameters as the mean of the parameters resulted in 
various iterations. The errors on these parameters are 
derived by determining the standard deviations. Fig. 
7 shows the histogram of the result obtained using the 
heuristic scan method for the orbital inclination, i in the 
contact mode with PHOEBE for the FO Hya. 

The light curve synthesis solutions from the present 
data indicate need of a third light. The presence of 
third light shows a tight correlation with the inclination 
angle and the mass ratio. It is remarked that correla- 
tion between mass ratio and third light still remains to 
be determined even for the best solutions available for 
the contact binaries showing total eclipses, where inner 
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eclipse contacts can be well defined and the mass ratio 
can be determined solely from light curves (Mochancki 
and Doughty, 1972). During light curve modeling of the 
FO Hya the mass ratio (q spec ) was fixed, therefore the 
third light highly correlated with the orbital inclination 
angle as shown in Fig. 7. In order to detect third light 
from photometric light curves with precision of a;0.01 
mag spectroscopic determination of the mass ratio is 
necessary (Pribulla and Rucinski, 2006). Precession of 
the orbital plane and an apsidal motion of the close bi- 
naries could be due to the third body orbiting a binary 
system. In eclipsing binaries, the change in orbital in- 
clinations and consequently, change in the photometric 
amplitude of the eclipses are due their precession of the 
orbital plane. 

The best fit parameters are given in the Table [4] 
The photometric elements such as T2, i, Qi, Q2, 
relative luminosity and spot parameters obtained in 
the present investigation agree well with the previous 
solutions given by Candy and Candy (1997) and Siwak 
et al. (2010). The fractional stellar radii r^2 obtained 
from the light curve modeling are normalised to the 
semi-major axis of the binary system, i.e., r = R/a, 
where a and R are the semi-major axis and component 
radii, respectively. The component radius for each star 
(in solar radius units) may be determined using the 
relation R12 = ri^a, where r is the geometrical mean 
of the r-pole, r-side and r-back radii for the FO Hya 
contact system. From the modeling of the light curves, 
the temperature difference between the components 
of FO Hya is found to be ~ 1787 K. Because of this 
large temperature difference the light curves of FO Hya 
show unequal depth. Such type of systems are known 
as B-type systems, which are in geometrical contact 
but not in thermal contact (e.g. Lucy and Wilson, 
1979). These binaries show an EB (J3 Lyrae)-type light 
curve. However, their orbital period falls in the range of 
classical W UMa systems (Kaluzny, 1986). These type 
of systems are defined as B-type contact systems which 
show large differences in depths of eclipses. These 
systems are yet appear to be in contact and exhibit 
P Lyrae-type light curves (Rucinski 1997). Some 
of them may be very close, semi-detached binaries 
mimicking contact systems and some may be in contact. 

We determine the absolute physical parameters of the 
components, based on the results of the light curve so- 
lution. The luminosity of the star is calculated using the 
temperature and radius of each of the individual com- 
ponents. Table|5]lists the results regarding the physical 
parameters of the individual components of the FO Hya 
system. 
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Figure 6: Geometric configurations of FO Hya generated by PHOEBE 
at phases (0) = 0.00, 0.25, 0.50 and 0.75, respectively. 



The eclipsing binaries are classified into contact, 
semi-detached and detached binaries according to fill 
factor (f) as < f < 1, f = and f < 0, respectively. 
The fill-out factor (f) is given by 



/ = 



Q,„ - Q 



1.2 



(12) 



where Q,„ and Cl out , and Q12 are the inner and outer 
Lagrangian surface potentials, and surface potentials of 
the star 1 and 2, respectively. In case of contact bi- 
naries surface potentials of primary and secondary are 
equal to the surface potential (O) of the common en- 
velope for the binary system (i.e. fii = Q.2 — Q). For 
semi-detached binaries with the primary and secondary 
components filling their Roche lobe, Qi = Q,„ and Q.2 = 
£2,„, respectively. For detached binaries Q\ .2 > £2,„. The 
fill-out factor for the FO Hya is obtained to be 0.68 and 
is in the range of contact binaries. 

The computed absolute parameters of FO Hya, as re- 
ported in Table [5] are used to estimate the evolution- 
ary status of the system by means of the mass-radius 
diagram shown in Fig. 8. In this plot, the solid line 
represents the zero age main sequence (ZAMS) from 
Schaifers and Voigt, (1982). The position of compo- 
nents on the mass-radius diagram are obviously outside 
the zero age main sequence and they are located left 
on the zero age main sequence (ZAMS). It is evident 
from the Fig. 8 that the radius of the primary compo- 
nent approaches to the value for zero age main sequence 
(ZAMS) stars, i.e. basically a ZAMS star (Awadalla and 
Hanna, 2005; Liu and Yang, 2000b). Fig. 8 shows the 
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Figure 7: Histograms of the results obtained using the Monte Carlo 
parameter scan for the parameters fit of inclination in the contact mode 
with PHOEBE for the FO Hya. 



Table 5: Absolute parameters obtained from the light curve modeling. 



Parameters 


Symbol 


Photometric 






elements 


Mass 






Primary 


(Mi/Ma) 


1.31+0.07 


Secondary 


(Ma/Me) 


0.31+0.11 


Radius 






Primary 


(Mean Radius/R Q ) 


1.62+0.03 


Secondary 


(Mean Radius/R Q ) 


0.91+0.02 


Luminosity 






Primary 


(Li/L ) 


5.65+0.21 


Secondary 


(L 2 /L ) 


0.55+0.03 



mass-radius relation of the secondary is different from 
that of ZAMS stars (Yang and Liu, 2001), i.e. the sec- 
ondary component is clearly far above the ZAMS than 
the location of primary, due to the fact that the mass 
transfer between the components will restructure the 
secondary and make it over-sized and over-luminous for 
its mass (Webbink, 2003; Li et al., 2008). 

7. Sources of Polarization 

The values of polarization and polarization position 
angle of FO Hya observed at three different epochs are 
given in Table [6] Most of the polarimetric observations 
of FO Hya have been taken during the phase from 
0.13-0.26. However, the polarimetric observations 
in B band have been taken during the phase of 0.64. 
The degrees of polarization are found well with in lcr 
during these phases. Thus the observed polarization 



0.5 1 1.5 2 2.5 3 

M/M ( , 

Figure 8: Location of the primary and secondary components of FO 
Hya on a mass-radius diagram. The continuous line shows the zero 
age main sequence. 



seems invariable. Since present observations are taken 
just before primary or secondary maxima so observed 
polarization may be invariable. The mean values of 
degree of polarization and polarization position angle 
were calculated to be 0.18+0.03, 0.15+0.03, 0.17+0.02 
and 0.15+0.02 per cent and 61+6, 115±6, 89+3 and 
98 ±5 degree in B, V, R and I bands, respectively. The 
degree of polarization in W UMa-type binaries were 
found upto 0.6% (e.g. Oshchepkov, 1978; Pirrola, 
1977). Oshchepkov (1978) found the variable polar- 
ization in W UMa type binaries and suggested that 
the photospheric scattering, the reflection effect and 
scattering in a gaseous envelope or stream could be 
the possible sources of variable polarization. However, 
Pirrola (1977) did not observe significant variability in 
the polarization and suggested that the polarization due 
to gas streams may be a transient phenomenon. More 
phase locked observations are needed to explain the 
observed polarization in FO Hya. 



Manset and Bastien (2000, 2001) have studied the 
polarimetric variations produced in a binary star by 
considering Thomson scattering and concluded that the 
average polarization produced depends on the grain's 
composition and size and on the wavelength of obser- 
vation and they find that the highest polarizations are 
produced by grains with sizes in the range ~ 0. 1 -0.2 yum. 
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Table 6: Observed BVRI polarization values for FO Hya. 



Date of 


Filter 


P) 


e 


JD 


Phase 


observations 




(%) 


O 


(2455000+) 




April 07, 2010 


B 


0.14±0.09 


62±11 


294.15901 


0.1752 




V 


0.17±0.02 


55±3 


294.16441 


0.1867 




R 


0.12±0.05 


149±6 


294.17191 


0.2027 




I 


0.11 ±0.04 


18±5 


294.16705 


0.1923 


April 23, 2010 


V 


0.12±0.08 


175±11 


310.16255 


0.2575 




R 


0.21±0.01 


28±1 


310.12981 


0.1878 






0.18±0.06 


178±7 


310.10323 


0.1312 


May 20, 2010 


B 


0.22±0.01 


60±1 


337.10598 


0.6382 



8. Summary 

We have carried out a detailed photometric analysis 
of the eclipsing binary FO Hya using the data obtained 
from B, V, R and I band CCD photometric observations. 
Making use of the data available in the literature our 
analysis is based on the data for 65 years. Our obser- 
vations yield the period to be 0.469556 ±0.000003 day, 
which is similar to that obtained previously. On the ba- 
sis of the long term data, we found that the period of FO 
Hya increases by 5.77 x 10~ 8 days yr~', probably caused 
by the mass transfer between the components at the rate 
of 1.67 xl0~ 8 M yr _1 . In the O-C curve analysis, it 
is also found that the period of FO Hya has varied in a 
sinusoidal way, superimposed on the long-term, upward 
parabolic variation. The period of the sinusoidal vari- 
ation is about 42.30 years with an amplitude of 0.004 
days. Further, the mass function and the mass of the 
third body are found to be 0.0002 M and 0. 10 M Q , re- 
spectively. The distance of FO Hya is determined to be 
446.68 pc. In the WD code analysis, we employed the 
mass ratio of 0.238±0.006 and found that the secondary 
component consist of dark spot on their surface. We 
have inferred the mass, radius and luminosity of pri- 
mary and secondary components to be 1.31 ±0.07 M 
and 0.31 ±0.11 M , 1.62±0.03 R Q and 0.91 ±0.02 Rq, 
and 5.65±0.21 L G and 0.55±0.03 L Q , respectively. We 
have also derived filling factor to be 0.68. In this analy- 
sis we have taken up, presumably for the first time, the 
polarimetric studies of FO Hya, and obtain the mean 
values of degree of polarization and polarization posi- 
tion angle to be 0.18±0.03, 0.15±0.03, 0.17±0.02 and 
0.15±0.02 per cent and 61+6, 115+6, 89±3 and 98±5 
degree in B, V, R and I bands, respectively. The level 
of optical polarization FO Hya could be due to the pres- 
ence of supplementary sources e.g. binarity and circum- 
stellar material (gas streams). 
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